International Atomic Energy Agency 
and 

United Nations Educational, Scientific and Cultural Organization 



INTERNATIONAL CENTRE 

FOR 

THEORETICAL PHYSICS 



Continuous Local Symmetry in Ising-type 

Models 

S. Randjbar-Daemi 
and 
J. Strathdee 



International Centre for Theoretical Physics 

January 1994 



Continuous Local Symmetry in Ising-type 

Models 



S. Randjbar-Daemi 

International Centre for Theoretical Physics, Trieste 34100, Italy 

and 

J. Strathdee 

International Centre for Theoretical Physics, Trieste 34100, Italy 

January 1994 

Abstract 

A class of generalized Ising models is examined with a view to extracting a low 
energy sector comprising Dirac fermions coupled to Yang-Mills vectors. The main 
feature of this approach is a set of gap equations, covariant with respect to one of 
the 4-dimensional crystallographic space groups. 
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Many authors have suggested that gauge fields may not be elementary In this 
note we would like to reconsider this possibility in the context of a discrete Ising-type 
model. We argue that it is feasible to set up a self-consistent scheme from which both 
gauge vectors and Dirac fermions emerge as effective dynamical variables in the continuum 
limit. 

Models of the Ising type would seem to provide the most primitive kind of structure 
on which one could hope to develop a semi-realistic low energy physics. We shall therefore 
consider an array of sites, A, on each of which there is an Ising variable S(A) = 0, 1. The 
partition function is given by the sum over configurations 

Z = J2exp f -f £ S(A)J(A, B)S(B) S(A)) (1) 

S \ Z A,B A J 

where the couplings, J(A, B) are to be specified. We have in mind a 4-dimensional lattice 
in which the couplings are not necessarily restricted to nearest neighbours. To restrict 
their form we shall assume that these couplings are invariant with respect to one of the 4- 
dimensional crystallographic space groups. In addition, we shall suppose that the sites are 
grouped into N- component cells such that the couplings are sensitive only to the cell and 
do not distinguish among the individual sites in a cell. This introduces a local permutation 
symmetry, Sjv, into the model. With these assumptions, the partition function flU) can 
be represented in the weak coupling, or high temprature, limit in the form of an integral 
over Grassmann variables, ipi(n), ipi(n), i — 1, ...,N, n e Z 4 , 

Z= ( (dipdip^j exp ( — i[j{n)ip{n) J(n — m)'0(m)'0(m) J (2) 

where the coupling function, J(n), is invariant under the action of the crystal point group, 
n — > gn. Equation (0) is the infinite temperature limit of a fermionic partition function. 
The external magnetic field ^ has been introduced in (|I|) in order to render the boundary 
conditions of the finite temperature fermionic path integral periodic. In the Grassmann 
version (fj) the local permutation symmetry is enlarged to local GL(N,C). 

Our purpose is to study the low energy properties of the system @. In particular, we 
want to know if there is a non- vanishing condensate < ip(n)ij)(m) > and if so, what are its 



symmetries. In order to obtain fermionic excitations carrying spin 1/2 it is essential that 
the ground state symmetry should involve an embedding of the point group of the lattice 
into the local group, GL(N,C), such that the Grassmann fields transform as spinors. It 
is for this reason that we assume the point group to belong to one of the crystal classes, 
i.e. be a discrete subgroup of 0(4). 

The largest 4-dimensional crystallographic point group, H, has 1152 elements ( in- 
cluding improper elements) || Q. Its spinor covering possesses a 4-dimensional spinor 
among its irreducible representations. We shall therefore assume that N is a multiple of 
4 and that the ground state symmetry is H x U(N/4) together with the lattice transla- 
tions. Under this group the ^-component if) decomposes into "Dirac" spinors of H, in 
the fundamental representation of U(N/4). 

To study the condensation problem we formulate self-consistent, Hartree-type equa- 
tions for the condensate amplitude 

c(n — m) =< if)(n)if)(m) > (3) 

and require the solution to be invariant with respect to H x U(N/4). This means, in 
particular the functions c(s) should satisfy 

c{gs) = a(g)c(s)a(g~ l ) (4) 

where g is an element of the point group and a(g) is its image in the local group, GL(N, C). 

A convenient way to obtain the self-consistency conditions is by introducing a bosonic 
auxiliary field, x(n, m), and replacing (0) by the equivalent form, 



I Uipdip) (d\) exp 



if)(n)x{n,m)ip(m) + 2j\n-m) 

(5) 

tr (x(n,m)x(m,n)) 



-In Det X En,m 2J(jl - m ) 



= J{dx) exp 

where the latter version is obtained by integrating the Grassmann variables. In the 
Hartree approximation, the condensate amplitude, c(n — m) =< x(n, m) >, is given by 
the saddle points of (||). These points are obtained by solving the "gap equations", 

c(s) = J(s)G(s) (6) 



1 The corresponding crystal lattice turns out to be the root lattice of the exceptional algebra F4 (3). 



where G is the fermion propagator. This propagator is to be expressed in terms of c(s), 

G(n - m) = J A* G(k) e^" m ) (7) 

where the integral extends over a cell of volume (27r) 4 , a Brillouin zone. The inverse of G 
is given by 

G-\k) = £c( S )e"^ (8) 

It has the same invariance properties as c(s) since the lattice sum in (§) is invariant under 
the action of H . Under this action the lattice decomposes into non-intersecting orbits. 
On each orbit J(s) is a constant and c(s) is expressible in terms of a finite number of 
parameters. These parameters are determined as functions of the couplings by solving 
the gap equations (Q) or, equivalently, by searching for extrema of the energy density 

We have verified that for the nearby orbits, s 2 = 2, 4 and 6, the invariance condition 
(|) implies 

c(s) = ci(s) + Ca^VTV* (9) 

where C\ and C2 are orbit invariants and 7^ are Dirac matrices. For more distant orbits 
there would be some proliferation of parameters. However, insofar as only the nearby 
orbits are involved, the general form of the fermion propagator is given by 

G-\k) = M(k)+i^ (10) 

where M and represent the invariant sums, 

M{k) =E s ci(s) e~* ks , Q^k) = ^E s C2(s)e- lks (11) 

It is known that the group H allows four algebraically independent invariants to 
be constructed from a single 4-vector, They are of order k 2 , fc 6 , k 8 and k 12 . Since k 2 
is also 0(4) invariant, it follows that G~ l (k) is 0(4) invariant in the neighbourhood of 
kn = if terms of order k 6 are neglected, i.e. 

G-\k) = Z^im + i jk) (12) 
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where Z^" 1 and m are linear in the condensate parameters c\ and c%. When the gap 
equations are solved they will be expressed in terms of the J' s. One can hope that the 
fermion mass, m, can be made small by tuning these couplings. 

The condensation phenomenon, which yielded propagating Dirac fermions, will also 
give rise to propagating bosons. These bosons are associated with the fluctuations of the 
"auxiliary" field, x(n, m). In particular, because of the local symmetry one expects to 
find Yang-Mills vectors by means of the ansatz, 



X{x,x 



s) = c(s) T ( exp i j:_ s a) (13) 



which should be valid in the continuum limit. In this limit, the vector transforms in 
the usual way, 

dA^x) = dj{x) -i[A^x),6{x)\ (14) 

where 9{x) belongs to the algebra of U(N/4), commutes with c(s), and is slowly varying 
on the lattice scale, \d^9\ << \6\. On substituting the ansatz flT3"| ) into the bosonic action 
(|), expanding in powers of A, dA,..., the leading terms must compose themselves into 
the standard Yang-Mills form. To identify the gauge coupling constant it is sufficient to 
examine the bilinear term 

\j -0ji E A%-k)A£(k)JiS(k) (15) 

where the sum over a = 1, 2, .., (N/4) 2 accounts for the local symmetry U(N/4). Invari- 
ance with respect to the gauge transformations (|i4|) implies 



Kl = -2 - Mv) + ... (16) 

near k = 0. To express the coupling constant , g 2 , in terms of the condensate parameters, 
Cx, c 2 , it would be necessary to expand the coefficients, 

A- 1 (fc) = E.,y Hf^s 



fj,u\ h '} — ^s,s' ^ 

iq(s+s') 



S -^tr (c(s)c(-s)) +J -0ptr (c(s)G(q)c(s')G(q - k)) e 



1 - e iks 



ks 



i— s v 



(17) 



around k = 0. 

The ansatz ([13]) can also be substituted into the fermionic part of the action defined 
by (|5p to obtain the coupling of the gauge field to the fermions. One finds the usual 
minimal coupling. To summarize, the low energy sector of the 4-dimensional generalized 
Ising model is describe by the effective action 

1 



S = d x 



a -i^i^ + Z a -V(/7 + m)V 



(18) 



where g 2 , Z^ 1 and m are ultimately to be expressed in terms of the lattice couplings, J(s), 
by solving the gap equations @ and substituting the resulting condensate parameters into 
(0) and (0). 

In this note we have attempted to minimize the number of assumptions about the 



Ising model coupling constants consistent with obtaining a low energy structure like (18 



In particular, we have not restricted the J' s to the nearest neighbours or next nearest , 
etc. Our reason for doing this is the hope that it will eventually be possible to treat the 
J' s as dynamical variables : quantities to be obtained by solving equations of motion. 
The 4-dimensional lattice structure should then emerge as a feature of the ground state 
symmetry. The space group would then have a status analogous to that of the Poincare 
group in general relativity. Indeed, such a theory would have to yield a long range tensorial 
interaction somewhat like gravitation. 
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